Techniques for Calculating two- loop Diagrams 



J. Fleischer H and 0. L. Vereti 



Fakultdt fur Physik, Universitdt Bielefeld, D- 336 15 Bielefeld, Germany. 

Abstract 

Methods developed by the Bielefeld-DESY-Dubna collaboration in 
recent years are: DIANA (Diagram ANAlyser), a program to produce 
"FORM input" for Feynman diagrams, starting from the Feynman rules; 
methods to calculate scalar diagrams: Taylor expansion in small mo- 
menta squared in connection with a mapping and the Pade method to 
sum the series. Recently program packages for the large mass expansion 
were written and applied to the Z — > bb decay. Reviews of these activities 
were presented in the proceedings of the Ustrori '97 and Rheinsberg '98 
conferences. Here we concentrate on recent developements in the large 
mass expansion, applied to the two-loop contribution of the Z — > bb de- 
cay in the mj, = approximation, taking into account higher order terms 
of the expansion in M^/mf. 

1 Introduction 

The calculation of diagrams with one non-zero external momentum squared 
(q 2 ) has wide applications in QED and QCD for both selfenergies and vertices. 
In these cases also only one non-zero mass enters the problem. In electroweak 
problems like Z — > bb one has mixing terms between electroweak and strong 
interactions and due to that different internal masses occur, so that the method 
of Taylor expansion is getting more difficult to apply. In this case, however, 
the top quark (m t ) plays a special role and it allows to make the expansion 
in the large mass. The method is not applicable to arbitrary high q 2 , but 
as has been demonstrated in El, for q 2 = m| this approach is still reliable. 
While in EJ only scalar diagrams have been considerd, here we investigate the 
full decay amplitude. It turns out that the obtained results are simpler for 
the full process than for scalar diagrams in the following sense: first of all, 
complicated functions like higher poly logarithms, which show up in the analytic 
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evaluation of scalar diagrams, cancel in the full amplitude; furthermore the 
convergence of the large mass expansion also turns out to be better for the full 
amplitude than for scalar diagrams. These observations are manifestation of 
gauge cancellations observed in gauge theories in general. Nevertheless their 
observation in this special form is surprising! 

Due to the fact that the method of mapping and Taylor expansion is quite 
useful and finds applications by other authors (see e.g. □), we give a short review 
here concerning the method and report on recent developements of calculating 
Taylor coefficients for two-loop diagrams. Then, in the second part, we turn 
to our main point namely the large mass expansion for the Z — > bb problem 
and the comparison with the work ofu. 

2 Expansion of three-point functions in terms of an external mo- 
mentum squared 

Taylor series expansions in terms of one external momentum squared, q 2 say, 
were considered for selfenergy diagrams incl, Pade approximants were intro- 
duced inO and in Ref. Eil it was demonstrated that this approach can be used 
to calculate Feynman diagrams on their cut by analytic continuation. In the 
case of a three-point function like Z — > bb in the mi, = limit we have for the 
external b— quark momenta p\ = p 2 = 0. The expansion of the scalar diagram 
then looks like 

oo 

c(pi, P 2) = J2 a n(pm) n , (i) 

with q 2 = {jpi +p 2 ) 2 - 

For the calculation xtf the Taylor coefficients in general various procedures 
have been proposed uau. These methods are well suited [far programming in 
terms of a formulae manipulating language like FORM Such programs, 
however, yield acceptable analytic results only in cases when not too many 
parameters (like masses) enter the problem. Otherwise numerical methods are 
needed U. 

In the case of only one non-zero mass and only one external momentum 
squared, indeed the case with the least nontrivial parameters, for many di- 
agrams analytic expressions for the Taylor coefficients can be obtained. For 
recent references seeE3. 

For the purpose of calculating Feynman diagrams in the kinematical do- 
main of interest it is necessary to calculate them from the Taylor series an their 
cut. This is performed by analytic continuation in terms of a mappingH 



Assume, the following Taylor expansion of a scalar diagram or a particular 
amplitude is given C(pi,p2, ■ • •) = J2m=o a my m = f{y) and the function on 
the r.h.s. has a cut for y > yo- 

The method of evaluation of the original series consists in a first step in 
a conformal mapping of the cut plane into the unit circle and secondly the 
reexpansion of the function under consideration into a power series w.r.t. the 
new conformal variable. We use 



1 - yi - y/yo 
i + v 71 - y/w ' 



(2) 
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Figure 1: Conformal mapping of the complex y- plane 



into the o;-plane. 



By this conformal transformation, the y-plane, cut from yo to +oo, is 
mapped into the unit circle (see Fig.0) and the cut itself is mapped on its 
boundary, the upper semicircle corresponding to the upper side of the cut. 
The origin goes into the point to = 0. 

After conformal transformation it is suggestive to improve the convergence 
of the new series w.r.t. u> by applying the Pade method I 13 !! 14 !. A convenient 
technique for the evaluation of Pade approximations is the e-algorithm of li-3 
which allows one to evaluate the Pade approximants recursively. 

Generally speaking, the precision of results with this mapping and Pade is 
of the order of 3-4 decimals with 30 Taylor coefficients for timelike q 2 values a 
factor of approximately 100 times the lowest threshold value. For lower q 2 (a 
few times the threshold value) the precision is of the order of 10 decimals in 
quite many cases. The precision worsens near second nonzero thresholds. 

As a final remark we mention that for diagrams with zero thresholds new 
techniques have been developed. In fact terms of the form ln m (<7 2 ) have to 



be factorized, where m is the number of zero thresholds ofpthfi diagram The 
factors in front are then expanded in terms of Taylor seriesli3'E3. 



3 Large Mass Expansion (LME) 

As mentioned above, for the evaluation of diagrams with several different 
masses, one of which being large (like the top mass rrit), we use the general 
method of asymptotic expansion in large masses £3. For a given scalar graph 
G the expansion in large mass is given by the formula 



F G (q,M,m,e) 



7 



FG/-y(Q,m,e) o T q ~< , m -< F 1 (q 1 , M, m 7 , e), (3) 



where 7's are subgraphs involved in the asymptotic expansion, G/7 denotes 
shrinking of 7 to a point; is the Feynman integral corresponding to 7; 
7g 7! m 7 is the Taylor operator expanding the integrand in small masses {to 7 } 
and external momenta {<z 7 } of the subgraph 7 ; o stands for the convolution 
of the subgraph expansion with the integrand Fqi^. The sum goes over all 
subgraphs 7 which (a) contain all lines with large masses, and (b) are one- 
particle irreducible w.r.t. light lines. 

For the Z — > bb decay we have q 2 = Af§ for the on-shell Z's. Fig. 2 shows 
diagrams with two different masses on virtual lines, one of which a top. W 
and Z are the gauge bosons with masses Mw and Mz, respectively; <f> is the 
charged would-be Goldstone boson (we use the Feynman gauge); t and b are 
the t- and b-quarks. Fig. 3 shows subdiagrams needed in the expansion (3). 



p2 



Pi 



P 2 



Pi 



g.7 



g.y 



g.7 b b b 

W, 



Z 

(p2-pl) 



Z 

(p2-pl) 




case 7 case 7.1 case 7.2 



Figure 2: Two-loop diagrams with two different masses on internal lines arising in Zbb. 
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Figure 3: The structure of the LME, see explanations in the text. 

Finally the LME of the above diagrams has the following general form 

1 N " / M 2 \ ' / n 2 \ j m m 2 

E (SO (h) E^\M^S)ln^ 

(4) 

where m is the highest degree of divergence (ultraviolet, infrared, collinear) 
in the various contributions to the LME (m < 3 in the cases considered). 
Myy/m 2 and q 2 /m 2 are considered as small parameters. Ai^ t k are in general 
complicated functions of the arguments, i.e. they may contain logarithms and 
higher poly logarithms. 

In contrary to the work ofS we see no inconveniences in directly applying 
the above method and did not use any "continued expansion" . 

In the following we present somepjof our results of the LME for the full 
two-loop 0(aa s ) contribution. As in 13, we are interested only in the virtual 
effect of the top quark, which renders the decay of the Z boson into bottom 
quarks different from the one into other down-type quarks. Therefore in the 
following the quantity T^_ d = T^^ bb — T ] ^_ rdd is considered in two loop order, 
in which expression the counterterm contributions cancel. The superscript W 
means that only diagrams with virtual W bosons are included. The other part 
with the Z exchange makesjio discrimination between b- and (i-quarks and is 
calculated in order aa s inliJ. Our result reads 
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where T° is the Born decay rate, y = M^/AM Wl Iq — In cos 8^, L t w = 



ln(mf/M^) and („ 
introduced 



£(n) the Riemann ^-function. The following integral is 



1 (-1)" / ln"(l - ty) 

2 n! J y/T—t 



dt. 



(6) 
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In the above final result enters only 1\ which has the expansion 
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(7) 

Inserting this expansion in (0) , we fully agree with the result of Harlander 
et.al. as far as they have presented their result. Our result is more compact, 
however, and it is interesting to observe that, while higher poly logarithms oc- 
cur in the scalar integrals!^, in the full decay amplitude these and the higher /„ 



(also, however, expressible in terms of poly logarithms) cancel. The remaining 
I\ , expanded above, is merely a logarithm: I\ = 2 — a log((a + l)/(a — 1)), a = 
a/I - l/y. Thus we observe that the final result is much simpler than inter- 
mediate results from scalar diagrams. Moreover, the convergence of the series 
in terms of large masses is much better than for the series obtained for scalar 
diagrams □, which is demonstrated below. 

Our numerical results are as follows: x\ and xi being the 1-loop and 2-loop 
results, r — M^y/ntf, the large mass expansion is given in the form 
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0.2360 - 1.0236r - 1.7492r 2 - 1.4477r 3 - 0.1758r 4 + 1.2997r 5 



x 2 - 



2.1005r 6 + 1.7986r 7 + 0.5692r 8 - 0.9260r 9 - 1.6273r 10 
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(9) 



0.7167 - 1.6940r - 4.0920r 2 - 4.7543r 3 - 2.9143r 4 + 0.5688r 5 



3.7747r 6 + 4.8688r 7 + 3.6932r 8 + 4.5196r 9 + 24.090r 1 



(10) 



The first terms in x\ , X2 correspond to the leading m 2 /M 2 ^. For m t — 175GeV, 
M w = 80.33GeV and M z = 91.187GeV we obtain 
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-0.5045-0.0704+ —(1.1556 + 0.1285) 



(11) 



In each of these terms the leading term and the corrections are given separately, 
te that the leading m 2 /M^ term in order 0(aa s ) was obtained earlier in 

The observation is that the series for the full amplitude converges like r n 
while for the scalar diagram the convergence was like (4r)™. Accordingly the 
term of order r 4 gives only an error of order 0.5% while for the scalar diagrams 
the errors where of the order several % EJ. 
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